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I. INTRODUCTION 


A. GENERAL 

The term biofouling refers to the settlement, 
attachment, and growth of marine organisms on surfaces 
that man puts into the ocean. This process has a profound 
impact on naval operations due to the fouling of ships 
hulls, rudders, salt water piping systems, sonar domes, 
and the fouling and biodeterioration of harbor or pier 
structures. This problem has been estimated to cost the 
Navy several hundred million dollars a year (Woods 
Hole,1952; Fisher et al,1975) due in part to increased fuel 
requirements caused by the greater frictional drag of 
fouled ships hulls, the increased repair or replacement 
cost of piping and machinery damaged by fouling organisms, 
and the need to expend funds to continually remove tne 
fouling organisms from vessels. 

The principal method used to combat the problem of 
biofouling on the hulls of ships has been the application 
of antifouling paints. In general, such paints contain any 
of several metallic compounds which, as they leach out of 
the paint matrix, are toxic to fouling organisms. 

The testing of antifouling paints prior to their 
release for general use involves the use of sampling assay 


techniques. In simple terms, a number of substrates are 


v5 








painted with the paint to be tested and are deployed in the 
sea for some arbitrary length of time. Following the 
retrieval of these substrates, they are compared with non- 
toxic control surfaces that have been exposed for the same 
length of time and the efficacy of the antifouling paint in 
preventing the settlement of fouling organisms is 
ascertained. 
1. Sampling Design 
The proper number of plates to deploy for the 
EectumagsOr antifouling paints, or for that matter, the 
study of fouling organisms in general, has always been 
somewhat arbitrary. This is because the proper number of 
plates to deploy to sample the fouling population is a 
direct function of the degree of variability within the 
population. Despite the great volume of fouling research 
that has been conducted, the study and quantification of 
this variability has only very recently been attempted. 
2. Previous Research on Fouling Community Variability 
Moist of EtEhe information dealing with the 
variability of biofouling communities has been collected 
within the past decade and is often somewhat contradictory. 
This contradiction is sometimes caused by the particular 
descriptors of the fouling community (percent cover, 


species counts, etc.) the researchers used in the analysis 


Ofeeres variabllity. 
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Research concerning the variability of small 
numbers of panels Suspended for only one month off the 
Florida coast (Mook,1976) showed very little variability in 
Species count. These results should be interpreted 
cautiously however due to the short time of immersion. 

Similar research conducted in North Carolina 
(Sutherland,1974) using a more extensive series of panels 
suggested that the development of the fouling community was 
extremely variable. This research supported the conclusion 
of earlier studies in California (Boyd,1972) which also 
found significant fouling community variability. 

Studies conducted in Hawaii (Schoener et al,1978) 
and in Massachusetts (Osman,1977) found that fouling 
community variability based on species counts was 
relatively low. This conclusion was echoed by a study which 
analyzed the variability of identical panels in terms of 
total percent cover, species count, and inter-panel 
Similarity indices (Schoener and Greene,1980). The results 
of the study indicated that approximately ten replicate 
panels were sufficient to resolve to a high degree of 
confidence, the mean value of these descriptors. 

As is evident, there is wide disparity between the 
various studies conducted to date on what the variability 
of the biofouling community at various locations truly is. 


Before meaningful antifouling paint test procedures can be 
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developed, proper sampling techniques based on the 
Milanmetetcacion Of fouling community variability must be 


devised. 


Be OBJECTIVE 

The primary objective of this thesis was to determine 
the variability of the biofouling community in Monterey 
Bay. Once this had been completed, the development of a 
appropriate sampling strategy based on this variability 
could be accomplished. This information was to be provided 
to the David Taylor Naval Ship Research and Development 
Center so that modifications to present antifouling test 


procedures could be undertaken as required. 
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Il. METHODS AND MATERIALS 


A. GENERAL 

The experiment consisted of deploying one-hundred (100) 
mild steel plates, each painted with one of four vinyl or 
epoxy based paints, in Monterey Harbor (Figure 1) for 
periods of up to 11 months and determining the variability 
of the resulting fouling communities that settled. Three of 
the paints used in the study contained antifouling 
compounds while the fourth ( the control surface ) did not. 
The fouling community structure on each plate was deter- 
mined destructively (the plates were not redeployed after 
Study) by microscopic analysis and the fouling population 


and makeup were determined. 


B. PLATES AND PAINTS 

The one-hundred plates used in this study were 
fabricated from low carbon, mild (structural) steel with 
dimensions 25.4cm x 30.5cm x .16cm. A small (.64cm) hole 
was drilled approximately 1.3 cm down from the midpoint of 
the top edge (one of the edges with the lesser dimension) 
of the plates to allow for the attachment of the suspending 
Pinemeore deo Woyment. Additionally, a small (.32cm x 2em) 
groove was milled through the plate approximately 6.5cm to 


One side of the drilled hole for the attachment of the 
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Figure 1. A Diagram of Monterey Bay Showing the Deployment 
Site at the Coast Guard Floating Dock. 
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identification tag. For consistent reference, the front of 
the plate was chosen as that side which would face the 
observer when the plate was held vertically with the dril- 
led hole at the top and the identification groove to its' 
fete (Figure 2). 

The four paints used in the experiment were: 

1. MICRON 22; a commercially available antifouling 
BommeecOomtalining Dbis{tributyltin) oxide and cuprous 
thiocyanate as the antifouling agents (Appendix A). 

2. Navy Standard Formula 121 Red Vinyl Antifouling 
Paeme, 8thne discontinued U.S. Navy antifouling paint 
containing cuprous oxide as the antifouling agent 
(Appendix B). 

3. Navy Standard Formula 170 Black Camoflage Vinyl 
MMmereouLing Paint; the currently used standard 
Peer OulLing paint of the U.S. Navy containing 
Pease TDuUcytin) oxide and tributyltin fluoride as the 
antifouling agents (Appendix C). 

4. Zynolyte Epoxy Rust Mate; a commercially 
available non-toxic corrosion resistant epoxy based paint 
used as the control surface (Appendix D). 

Allplates were sandblasted then primed and painted in 
accordance with label directions. 
1. Priming Procedure 
The priming procedure consisted of first applying 
one coat of Navy Standard Formula 117 'Green Wash Primer' 
to both sides of all plates with a 5cm latex rubber brush 
and allowing this to dry for twenty-four hours. This was 
followed by the application of two coats of Navy Standard 
Formula 119 'Red Lead' to both sides of all plates using a 


fscemenylon paimesroller. The first coat of Formula 119 was 
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Figure 2. The Front Side of One of the Experimental 
Plates. 
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allowed to dry for forty-eight hours before the second coat 
was applied, and an additional forty-eight hours drying 
time was allowed prior to painting with the vinyl or epoxy 
paints. 
2. Painting Procedure 

All painting was accomplished using 7.6cm nylon 
paint rollers. The plates were divided into four groups 
(A,B,C, and D) of twenty-five plates each and labelled and 
painted as follows: 


1. Plates Al through A25 were painted on both sides 
with one coat of MICRON 22 antifouling paint. 


2. Plates Bl through B25 were painted on both sides 
with one coat of Navy Standard Formula 121 Antifouling 
Paint. 

3. Plates Cl through C25 were painted on both sides 
with one coat of Navy Standard Formula 170 Antifouling 
paint. 


4. Plates Dl through D25 were painted on both sides 
with one coat of Zynolyte Epoxy Rust Mate. 


The paint on all plates was then allowed to cure 
for 96 hours prior to deployment. The plates were labelled 
sequentially in each group (Al,A2,etc) by affixing, through 
the identification groove, an embossed DYMO tape label to 
each plate using monofilament nylon fishing line. This 
method of labelling was chosen to prevent the catalytic 
corrosion problems attendant with standard bronze or copper 


tags in contact with the steel plates. 


23 








C. DEPLOYMENT PROCEDURE 

The plates were randomly divided into twelve groups of 
eight plates consisting of two plates from each of the four 
Paint groups. To each plate a one meter length of .32 cm 
diameter stainless steel cable was then affixed by passing 
the cable through the drilled hole and forming a loop which 
was closed using Nico-Press crimp fittings (Figure 3). A 
Similar loop was formed at the distal end of the cable to 
facilitate attachment at the deployment site. 

On 22 and 23 May 1982, the plates were suspended (by 
groups) beneath the service access covers that extend the 
length of the floating dock at the Coast Guard Station 
Monterey. Each plate was individually deployed by attaching 
the distal loop of the cable to 10d nails driven into the 
dock and allowing the plate to hang vertically in the water 
(Figure 4). The plates were separated by a minimum 
horizontal distance of 40cm with the tops of the plates 
approximately one-half meter beneath the surface of the 
water. Since the dock rose and fell with the tide, the 
depth of immersion remained constant. Water depths below 
the plates (at MLLW) ranged from approximately three meters 
depth at the shallower end of the dock to more than ten 
meters depth at the seaward end. After being submerged for 
one month, inspection of the plates revealed that 


significant galvanic corrosion had occured at the junction 
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Figure 3. A Drawing Showing the Method of Attachment of 
the Suspending Cable and Identification Tag. 
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Figure 4. A Perspective View in Cross Section Showing the 
Deployment of the Plates at the Coast Guard Dock 
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between the mild steel plates and the stainless steel 
cable. As a result, on 23 June 1982, all stainless steel 
cables were removed and replaced with .95 cm diameter nylon 
line. This necessitated the exposure of each plate to the 
atmosphere for approximately thirty seconds during the 
replacement operation but , since the plates remained 
moist, it was felt no harm was done to the fouling 


organisms that had settled. 


D. FOULING COMMUNITY CENSUS AND IDENTIFICATION 

Bach month following the initial deployment, one of the 
groups of eight plates was randomly selected for retrieval 
and study. 

1. Sampling Procedure 

For each group of plates, a sampling grid of one- 

hundred (100) uniformly distributed random points was 
generated and graphically plotted by computer on a 25.4 cm 
x 30.5 cm output sheet. These points were then transferred 
manually to a clear plexiglass cover. The fouling 
communities on the plates were then systematically analyzed 
by setting the plates horizontally in a shallow container 
filled with seawater and positioning the plexiglass cover 
over the top coincident with the edges of the plate. 
Animals beneath plotted points were then censused and 
identified through the use of a stereo microscope. Since 


the plates had been suspended vertically in the water, none 
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of the sedimentation problems associated with horizontal 
deployment strategies developed. Therefore both sides of 
the plates were analyzed and counted as separate 
substrates. 
2. Identification 

To eliminate the necessity for compound microscope 
identification of settled organisms and to negate the ef- 
fects of neuston contamination of census results, only 
sessile, attached organisms greater than .5mm in size were 
counted and identified. Identification was accomplished 
through the use of available keys and literature; 
(Osburn,1952: Knight Jones,1979; Haderlie,1974; Morris 
Peed; Frrazier,193/; Smith and Carlton,1975). In all, 
some thirty-two taxa were identified aS major space occu- 


piers during the study (see Appendix H for species list). 
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A. EXPERIMENTAL DATA 

The data obtained from the census and identification of 
the organisms on the plates retrieved each month was used 
Bemgqenmerate —the Statistics that described the fouling 
populations. The two main descriptive quantities used to 
assess the variability of the fouling community were the 
total percent fouling cover on each plate and the 
Similarity of the fouling organizations between plates. 

1. Percent Cover 

The initial estimate of the total percent fouling 

cover on each plate was calculated by dividing the number 
of points from the census that had organisms beneath them 
(this value was termed S,) by the total number of points 
censused per plate (Ny = 100 for all plates). Since this 
method was obviously subject to some unknown degree of 
uncertainty, revised estimates for the total percent 
fouling cover were made using statistical techniques 
described below. 

ZOE Micarity 

The similarity between the fouling communities on 

the plates retreived each month was determined by the 
erplecimecemon Of Che Bray-Curtis Similarity index, I. 


(Whitaker,1952). This index is defined as: 
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ee = }min(a,b) 
where a and b are the fractional species proportions 
Peescmt om plates A and B. These fractional species 
proportions were determined from the initial census data by 
dividing the total number of instances a particular species 
was counted during the plate census by 100. In this study, 
empty points not occupied by any organism were treated as a 
Separate species. As an example, suppose the following data 


were collected: 


PrOpOrt1On on PEODOrtiloneon 
Species Plate A (=a) Plate B (=b) min(a,b) 
Empty ow ~24 24 
Species #1 0 lal: 0 
Species #2 46 aa 42 
Species #3 24 23 20 
then the Similarity Index = I, = >) min(a,b) = .89 


The Similarity index is thus a measure of the degree of 
Variability between the fouling communities on the two 
plates in terms of both the total percent cover and the 
species composition. 

The main purpose of calculating the similarity 
index between plates was to examine whether any of the 
antifouling coated plates displayed more variability than 
did the non-toxic control surfacés. This was of particular 
interest since a previous thesis (Kelley,1981) which dealt 
with marine microfoulers in Monterey Bay had shown that 


Surfaces coated with an organo-metallic antifouling paint 
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served as attractants to the organisms and, hence, were 
more heavily fouled, showed greater species diversity , and 
were more variable in terms of their fouling communities 
than were the control surfaces. Provided that the 
macrofouling community investigated in this thesis did not 
behave similarly, that is if the fouling communities on the 
non-toxic control surfaces were consistently more variable 
than were those on the antifouling coated surfaces, then 
any sampling strategy which could discern to an acceptable 
degree of error the amount of variability of the 
communities on the control surfaces would be able to do at 
least as well concerning the antifouling coated surfaces. 
The thrust of this thesis was to first investigate whether 
or not the control surfaces for each month exhibited more 
variability than did the antifouling coated surfaces. 
Provided this criteria was met, the next step was to devise 
an appropriate sampling strategy for the control surfaces 
using sophisticated statistical modelling to extrapolate 
the data obtained from the four control surfaces retrieved 


each month to any desired number of simulated samples. 


B. STATISTICAL MODELLING 

Once the initial percent fouling coverage estimates for 
the four non-toxic control surfaces retreived each month 
were determined, the data had to be manipulated to obtain 


better estimates of the variability of the superpopulation 
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of marine fouling organisms the plates were assumed to 
be sampling. Since there were only four substrates examined 
per month, a model capable of extending the experimental 
data to simulated samples of any size was required. The 
model also had to be capable of allowing and quantifying 
inter-plate variability within the simulated sample set. 
The techniques used to successfully meet these requirements 
have only recently been developed and this thesis is 
apparently the first incorporation of these techniques into 
fouling research. 
1. Model Alternatives 

One of the methods that has been used in the past 
(Schoener and Greene, 1980) to examine the degree of 
variability of biofouling coverage values has been to 
assume that the fractional coverage estimates have a normal 
distribution. Using this simple model, the upper and lower 
95% confidence limits are calculated about the experimental 
mean percent fouling cover by use of the formula: 

the 95% Confidence Limits = ¥ [+/-] 1.96* @ /YN- 
Wiewemx is the mean percent fouling cover, @ is the 
experimental standard deviation, N is the number of samples 
examined , and the formula is derived from the standard 
normal distribution. By assuming that the mean (x) and 
Standard deviation (@) do not vary with the number of 


plates examined, one varies N in the formula to determine 
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the effect that the number of plates has on the approach of 
the upper and lower confidence limits to the experimental 
mean. It must be pointed out however that this model has 
serious drawbacks. The first problem with this approach is 
the assumption that fractional fouling coverage values 
which will always lie between zero and one can be described 
by a normal distribution that is unbounded in range. While 
this assumption permits the calculation of various statis- 
tical parameters uSing well known analytical formulas, it 
1s obviously weak theoretically. Secondly, this model 
requires that the mean and standard deviation are stable 
with increaSing sample size. This requirement has the 
effect of forbidding inter-plate variability. Not only is 
such a restriction biologically untenable, it reduces the 
formula for the calculation of the confidence limits to : 

the 95% Confidence Limits =(+/-] constant/N 
Note that the above equation will result in a curve similar 
to [+/-] 1/JN- when plotted. Since the value of 1/ JN 
decreases by nearly 70% as N goes from one to ten, use of 
this method will always show that approximately ten plates 
are sufficient to resolve the variabilityof the mean 
percent cover. This result however is merely an artifact of 
the simplistic model used in its calculation. 

Another possible method that could be used to 


estimate the degree of variability of the biofouling 
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community would be to invent a predictive model for fouling 
populations. Unfortunately, such a model is considerably 
beyond our abilities at the present time. Such a model 
would require advective/diffusive models accurate spatially 
to microscale resolution . It would also require the 
ability to parameterize the entire range of biological 
forcing functions which include predation, nutrient supply, 
Susceptibility to environmental fluctuations, behaviorism 
of the organisms involved and planktonic larva survival 
rates to name just a few. Obviously, such sophistication in 
a model is not likely in the forseeable future. 

Since a stochastic model which assumed the fouling 
coverage values had a normal distribution was insufficient 
due to its’ inherent restrictions which forbid interplate 
Variability, and since a predictive model was unattainable, 
this study used a probabilistic model to explore the 
varlability of the fouling community. In such a model, the 
variability of the initial experimental data is determined 
and the statistical descriptors of that varibility are used 
as inputs for computer simulations of the experiment. 

2. Procedure 

Batewe plate (t) was assumed to have a 

random,independent proportion of fouling, Pe PB was assumed 


tombe of the £f£erm: e 





where €é, was assumed to be normally distributed with unknown 
mean (uw) and variance (a°) which were independent of the 
other plates censused at the same time. Note that the above 
equation could be written as: 
€ = os 
a In(P, 7 Pe )) 

Bach month, the €, value For wecach, Of Ehe folr umtresnted 
plates was determined as were the mean and variance of the 
four €, values. Revised estimates for these parameters were 
then determined using the Method of Maximum Likelihood (see 
Appendix E for mathematical development). 

Using as inputs the maximum likelihood values for 

2 

the mean (i) and variance (&@) of the monthly epsilon 
values, bootstrap computer simulations (Efron,1979) of the 
experiment were conducted. The bootstrap method was used to 
Simulate 200 groups of various fixed numbers of plates (see 
Appendix F for a complete discussion). 

Using the results from the bootstrap computer 
Simulations, the following statistics were calculated: 

1. The expected value of the mean percent fouling 
coverage. This was determined by calculating the mean of 
the 200 simulated group means for each of the various 
numbers of plates per group simulated. 

2. The upper 95% confidence limit about the expected 
value of the mean percent fouling cover for each of the 
various numbers of plates per group simulated. This was 
estimated by ordering the mean percent foulng coverage 
values for each of the 200 simulated groups of plates in 


ascending order and then finding the 95% quantile of the 
expected mean percent fouling coverage values. 
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3. The percent fouling cover that one can say 90% of 
the simulated plates will have fouling coverage less than 
Or equal to (with 95% confidence). This was done by 
Finding the 95% quantile (as described above) of the 
parameter NEWCVR(I) where: 


B B 


NEWCVR(I) =e /(l+e ) 
and B = f(I) + 1.285* G (TI) 


Note that it (1) and @(I) are the bootstrap simulations 
of the mean and standard deviation of for each of the 
200 group simulations. The value 1.285 is the 903 
quantile for the standard normal distribution. 


The purpose of calculating the parameter NEWCVR was 
to estimate how many plates it would take to be 95% 
confident of capturing 90% of the variance of the 
biofouling population. The confidence and variance values 
chosen were arbitrary, and the model can be readily 
modified to estimate the number of plates required to 
ascertain any degree of variance to any desired 
confidence. 


4. The standard deviation about the mean percent 
fouling coverage that 90% of the simulated plates would 
have was also estimated . This was done by finding the 
mean and standard deviation of the 200 NEWCVR(I) values 
for each of the various numbers of plateS per group 


Simulated each month. 
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[V. RESvutsS 


A. GENERAL 
Experimental data were collected monthly and the 
initial percent cover and similarity values computed (and 
the coverage estimates tabulated) for the second through 
the eleventh months of the experiment (Appendix G). 
Following the procedure discussed in the statistics 
chapter, these initial estimates were used as starting 
values for Bootstrap/Maximum Likelihood computer 
Simulations of the experiment. The computer simulation 
results from month 2 are described in some detail. The 
Simulation uSing the data from months 3 through 11 followed 
the same procedure and are each presented in a brief 
synopsis. 
1. Method Verification 
The procedure described in this thesis to develop 
an appropriate sampling strategy for antifouling paint test 
purposes was based on the assumption that any strategy 
which could ascertain the variance of the biofouling 
community on non-toxic control surfaces would be able to do 
at least as well regarding antifouling coated surfaces. 
This assumption would be correct provided the fouling com- 
munities on non-toxic control surfaces were more variable 


than were those on the antifouling coated surfaces. By 
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analyzing the monthly mean similarity values for the 
fouling communities on these two types of surfaces (Figure 
Sy lt is clear that the antifouling coated surfaces dis- 
played consistently less variability in their fouling 
PaedeciuEe than did the control surfaces. Therefore, this 
assumption is considered to be quite strong. 

A second major assumption used in the development 
of this procedure was that the monthly epsilon (€,) values 
Wsedeinm the calculation of the proportion of fouling oar 
were independent random values from a normal distribution 
with mean (uu) and variance ( Gan 

If the epsilon values were in fact distributed 
normally, then dividing the individual epsilon values for 
each month by the monthly standard deviation of the epsilon 
values and subtracting the mean monthly value of epsilon 
from this result, would transform the epsilon distribution 
into the standard normal distribution. By plotting the (N) 
transformed epsilon values ranked in ascending order 
against the theoretical order statistic obtained by finding 
the inverse function of the standard normal distribution 
for (j/N+l) as j goes from 1 to N, a diagnostic plot of the 
model was obtained. If the model was perfect, there would 
be a one-to-one correspondence between the transformed 
epsilon values and the order statistic. Plotting the forty 


transformed epsilon values obtained in this study against 
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Figure 5. Chart Showing the Monthly Mean Similarity Values 
for the Control Surfaces (Solid) and the Antifouling Coated 
Surfaces (Dashed). 
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the order statistic (Figure 6) shows that the assumption of 
a normal distribution for the variable epsilon is quite 
good. The least squares best fit line for the data had a 
ceeope Of 1.07 (vice 1.00 for the theoretical perfect 
correspondence) and the correlation coefficient between the 
transformed epsilon values and the order statistic was 
O98. 

Since the initial experimental data for percent 
fouling cover and the Bootstrap/Maximum Likelihood model 
were not independent but rather were coupled by the esti- 
mates of the mean and variance of the monthly epsilon 
values, one would reasonably expect the final mean percent 
fouling cover predicted by the Bootstrap simulations to 
approximate the actual data value if the model behaved 
reasonably. A chart comparing the arithmetic mean of the 
monthly values for percent cover obtained experimentally to 
the expected value of the mean percent cover predicted by 
the bootstrap simulations (Figure 7), shows that the model 
agrees quite well with the experimental data. In those 
instances where there was a Significant difference between 
the mean percent cover obtained from the data and that 
predicted by the bootstrap model, analysis of the actual 
data suggested that simply finding the arithmetic mean of 
the four monthly data percent coverage estimates was per- 


haps too sensitive to unusually high or low values. 
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Figure 6. A Diagnostic Plot of the Model With the Ranked, 
Normalized Values for the Individual Epsilon Values Plotted 
on the Ordinate (Abbreviated as Z) and the Theoretical 
Order Statistic Plotted of the Abscissa. The Dashed Line 
Indicates Perfect Correspondence and the Dotted Line is the 
Least Squares Best Fit for the Data. 
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Figure 7. Chart Showing the Arithmetic Mean of the Percent 


Fouling Cover From the Experimental Data (Solid) and the 


Bootstrap Simulated Mean Percent Cover (Dashed). 
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In summary, it is felt that the assumptions made in 
the formulation of this model are quite good and that the 
model gives very reasonable estimates of the most likely 
value for the percent fouling cover. 

2. Explanation of Figures 

The similarity indices between the plates retrieved 
each month were plotted on two separate graphs for the 
control surfaces and the antifouling coated surfaces. The 
range of the similarity values from zero to one was’ plot- 
ted on the ordinate and the plate designations were listed 
on the abscissa. The 'F' and 'B' that follow the plate 
group number refer respectively to the front and back of 
the plate. To obtain the similarity value between any two 
plates, find the horizontal line above the abscissa with 
arrowheads which terminate at points above the desired 
plate designations, then proceed horizontally to the 
ordinate to find the similarity index value. 

The expected value of the computer simulated mean 
percent fouling cover and the upper 95% confidence limit 
about the mean as functions of increasing sample size were 
plotted for each month. The expected value of the mean 
percent cover was dispayed as a solid line and the upper 
95% confidence limit was dashed. 

Graphs were also prepared showing the results of 


the monthly computer simulations used to determine the 95% 
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confidence limit of the percent fouling coverage that 903 
of the simulated plates will have fouling coverage less 


than or equal to. The ordinate of these figures was gra- 


duated in percent cover from zero to one hundred percent. 


The abscissa was labelled with the number of plates per 


group used in the simulations. The 95% confidence limit for 
the percent cover of 90% of the plates was then plotted as 
a function of increasing sample size. 

The same units for the ordinate and abscissa were 
used in the figures showing the standard deviation about 
the mean percent cover that 90% of the plates would have as 
Functions of increasing sample size. For these figures, 
however, the percent cover label on the ordinate refers to 
the deviation about the mean percent cover and not the 


total percent fouling cover these plates will have. 


B. RESULTS FROM MONTH 2 
1. Experimental Data 
The fouling community on the control surfaces after 
two monhs immersion was dominated by the hydroid Obelia 
spp. Ectoproct colonies each consisting of several dozen 


zooids had also settled by this time. These initial 


bryozoans were primarily Hippothoa hyalina and Celloporaria 


present.Watersipora is indigenous to the Galapagos Islands 
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and has not been reported in the literature farther north 
than the Gulf of California (Osburn,1952). It has not been 
noted in fouling studies conducted in Monterey Bay over the 
last twenty years and its appearance on the plates in this 
Study 1S probably due to the abnormally warm coastal waters 
caused by the recent El Nino event. 

The Similarity values for the control surfaces 
ranged from .74 to .92 with a mean of .83 (Figure 8A). The 
lack of any settlement on the antifouling coated surfaces 
resulted in similarity values of 1.0 for all of those 
plates (Figure 8B). 

2. Computer Simulations 

The bootstrap computer simulation of the expected 
mean percent fouling cover (Figure 9) resulted ina final 
iterated estimate of 44.4% for the simulated mean percent 
cover (vice 44.5% for the data arithmetic mean). Note that 
the dashed line indicating the upper 95% confidence limit 
of the mean percent cover does show some inverse 
relationship toN (the number of samples) but does not show 
the precipitous approach to the mean as N goes from one to 
ten that was predicted by the model which assumed a normal 
distribution for the fractional coverage values. 

Perhaps the most striking result of the computer 
Simulations was that of the effect of increasing the 


number of simulated plates sampled on the 95% confidence 
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Figure 8. Similarity Graphs for (A) Non-treated Control 
Surfaces and (B) Anti-fouling Coated Surfaces for Month 2. 
Dashed Lines Indicate Mean Similarity Values. 
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Figure 9. Computer Simulations Using Data From Month z 
Showing the Expected Value of the Mean Percent Fouling 
Cover (the Mean of the 200 Individual Group Simulation 
Percent Fouling Covers) as a Solid Line and the 95% 
Quantile (Dashed) of the Expected Mean Percent Fouling 
Cover as a Function of Increasing Sample Size. 
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limit of the percent fouling cover that 90% of the plates 
would have (Figure 10A). Note that there is very little 
dependence on N on the percent cover that one can say 90% 
of the plates will have fouling cover less than or equal 
Pemmetmais iS in spite of the fact that the standard 
deviation that these 90% of the plates will have about the 
mean percent cover does display a strong inverse dependence 
on N (Figure 10B). While these two results might seem 
dichotomous, they are not. The standard deviation about the 
mean does decrease with increasing sample size just as one 
would expect from the Central Limit Theorem. The fact that 
the 95% confidence limit of the fouling cover of 90% of the 
plates does not behave similarly is simply a result of the 
fact that by allowing interplate variability, one no longer 


constrains the parameter to have a 1l/ N dependence. 


C. RESULTS FROM MONTH 3 
1. Experimental Data 

The fouling community structure on the control 
surfaces after three months immersion remained dominated by 
the hydroid Obelia spp. The four species of spirorbid 
worms that live in Monterey Bay ,particularly Circeis 
armoricana, were also present. The protozoan Folliculina 
sp.was present in large numbers on three of the plates. The 
bryozoans were represented by four species with 


Celloporaria brunnea dominating this group. 
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Figure 10. Computer Simulations Using Data From Month 2 
Showing : (A) the Percent Fouling Coverage That One Can Say 
90% of the Samples Will Have Fouling Coverage Less Than or 
Equal to { With 95% Confidence); and (B) the Standard 
Deviation of Percent Coverage (About the Mean) That These 
samples Will Have As a Function of Increasing Sample Size. 
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The Similarity values ranged from .57 to .86 with a 
Meal OL ./4 for the control surfaces (Figure IIA). The 
mean Similarity for the antifouling coated surfaces 
Bemained at 1.0 (Figure 11B) since ‘Wo settlement of 
organisms had occured. 

2. Computer Simulations 

The computer simulations of the expected value of 
the mean percent fouling cover (Figure 12) resulted in a 
final iterated value of 49.5% vice a 49.3% arithmetic mean 
of the initial data. The upper 95% confidence limit about 
the mean decreased by nearly 8% as the number of simulated 
plates per group increased from 2 to 15 and then only 
decreased an additional 3% as the number of plates per 
group wasS increased to 80. 

The computer simulations of the 95% confidence 
limit of the fouling cover that 90% of the plates would be 
Fouled to an equal or lesser extent showed the value 
remained nearly constant at approximately 70% (Figure 
13A).The standard deviation of the percent cover on these 
plates about the mean percent cover showed little relative 
decrease past the simulated analysis of 20 plates per group 


(Figure 13B). 
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Figure ll. Similarity Graphs for (A) Non-treated Control 
Surfaces and (B) Anti-fouling Coated Surfaces for Month a 
Dashed Lines Indicate Mean Similarity Values. 
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Figure 12. Computer Simulations Using Data From Month 3 
Showing the Expected Value of the Mean Percent Fouling 
Cover (the Mean of the 200 Individual Group Simulation 
Percent Fouling Covers) as a Solid Line and the 958% 
Quantile (Dashed) of the Expected Mean Percent Fouling 
Cover as a Function of Increasing Sample Size. 
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Figure 13. Computer Simulations Using Data From Month 3 
Showing : (A) the Percent Fouling Coverage That One Can Say 
90% of the Samples Will Have Fouling Coverage Less Than or 
Equal to ( With 95% Confidence); and (B) the Standard 
Deviation of Percent Coverage (About the Mean) That These 
Samples Will Have As a Function of Increasing Sample Size. 
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D. RESULTS FROM MONTH 4 
l. Experimental Data 

ime LOuling community On the control surfaces after 
four months immersion was dominated by the bryozoan 
Watersipora cuculatta. Thehydroid Obelia spp. had been 
overgrown to a large extent by W. cuculatta and as a result 
its relative contribution to the percent cover was 
diminished considerably. In addition to W. cuculatta, there 
were Six other species of either upright or encrusting 
bryozoans present. Various spirorbid and serpulid worms 
were present in limited numbers. 

The percent coverage estimates varied greatly from 
plate to plate with a maximum of 84% and a minimum estimate 
on one plate of 2% cover. The wide variability was perhaps 
partially caused by the ascendency of the bryozoans as the 
dominant organism but the reason for the nearly total lack 
of fouling on two of the surfaces is unknown. 

The wide range in the percent coverage estimates 
was mirrored in the variability of the similarity indices 
between the control surfaces. The values ranged from .16 to 
.98 (Figure 14A). The antifouling coated plates still had 
no fouling so again the mean Similarity was 1.0 for those 


plates. 
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Figure 14. Similarity Graphs for (A) Non-treated Control 


Surfaces and (B) Anti~fouling Coated Surfaces for Month 4. 
Dashed Lines Indicate Mean Similarity Values. 
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2. Computer Simulations 

The final iterated value for the simulated mean 
percent fouling cover was 31% (Figure 15). The upper 95% 
confidence limit about this mean was quite large and 
decreased by nearly 40% as the number of simulated plates 
per group went from two to twenty (Figure 15). 

The wide variability of the percent coverage 
estimates caused the 95% confidence about the fouling 
coverage estimate of 90% of the plates to be quite high 
(Figure 16A). As can be seen, the only statement that can 
be made about the fouling coverage that 90% of the plates 
willhave is that the coverage will be something less than 
93% (Figure 16A) even though the estimated mean percent 
coverage was 31%. 

The standard deviation that these 90% of the plates 
will have about the mean was also quite large (Figure 16B) 
and showed significant decline out to approximately the 


thirty plates per group simulation point. 


E. RESULTS FROM MONTH 5 


1. Experimental Data 
The bryozoan dominance of the fouling community was 
firmly established by the fifth month of immersion. 
Watersipora cuculatta was the primary fouler on three of 


the four control surfaces and there were an additional 
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Figure 15. Computer Simulations Using Data From Month 4 
Showing the Expected Value of the Mean Percent Fouling 
Cover (the Mean of the 200 Individual Group Simulation 
Percent Fouling Covers) as a Solid Line and the 95% 
Quantile (Dashed) of the Expected Mean Percent Fouling 
Cover as a Function of IncreaSing Sample Size. 


ay 








73.0 100.0 


90.0 


(A) 


PERCENT COVER 


25.0 





0.0 


0.0 5.0°10.0 15.0 20.0 25.0 30.0 35.0 40.0 45.0 50.0 55.0 60.0 65.0 70.0 75.0 80.0 
NO. OF PLATES 


CO 


PERCENT COVER 
15,0) 920.0" 25.0 3850 35.07 40.0 


10.0 


0.0 5.0 


0.05.0 10.0 15.0 20.0 25.0 30.0 35.0 40.0 45.0 50.0 55.0 60.0 65.0 70.C 75.0 80.0 
NSE 1elintes: 


Figure 16. Computer Simulations Using Data From Month 4 
Showing : (A) the Percent Fouling Coverage That One Can say 
90% of the Samples Will Have Fouling Coverage Less Than or 
Equal to ( With 95% Confidence); and (B) the Standard 
Deviation of Percent Coverage (About the Mean) That These 
Samples Will Have As a Function of Increasing Sample Size. 
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eight other species of bryozoans identified during the 
census. 

maece was a fairly large range of initial percent 
cover estimates for the control surfaces with a minimum of 
13% and a maximum of 96%. The mean percent coverage 
estimate for these plates was 63% and the standard 
deviation was 31%. 

Pie ssimtlarity imdices for the control surtace 
plates ranged from .17 to .67 with a mean value of .45 
(Figure 17A). The similarityvalues For the 
antifouling coated surfaces ranged from .85 to 1.0 (Figure 
17B). This variability was caused by the settlement of 
Opelia spp. and the’spirorbid worm Circeis armoricana on 
several of the antifouling coated plates. 

2. Computer Simulations 

The final iterated estimate for the mean percent 
cover was 69% (Figure 18) vice the arithemetic mean of 63% 
for the initial data. The 95% coinfidence limit about this 
mean showed an appreciable decrease out to approximately 
twenty plates per group simulated. 

The population variability again resulted in the 
95% confidence limit of the fouling cover of 90% of the 
plates having virtually no dependence of the number of 
Samples analyzed (Figure 19A). The standard deviation about 


the mean that these 90% of the plates would have showed 
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Figure 17. Similarity Graphs for (A) Non-treated Control 


Surfaces and (B) Anti-fouling Coated Surfaces for Month 5. 
Dashed Lines Indicate Mean Similarity Values. 
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Figure 18. Computer Simulations Using Data From Month 5 
Showing the Expected Value of the Mean Percent Fouling 
Cover (the Mean of the 200 Individual Group Simulation 
Percent Fouling Covers) as a Solid Line and the 953 
Quantile (Dashed) of the Expected Mean Percent Fouling 
Cover as a Function of Increasing Sample Size. 
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Figure 19. Computer Simulations Using Data From Month 5 
Showing : (A) the Percent Fouling Coverage That One Can Say 
90% of the Samples Will Have Fouling Coverage Less Than or 
Equal to ( With 95% Confidence); and (B) the Standard 
Deviation of Percent Coverage (About the Mean) That These 
Samples Will Have As a Function of Increasing Sample Size. 
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little decrease past approximately 20 plates per group 


Simulated (Figure 198B). 


F. RESULTS FROM MONTH 6 

l. Experimental Data 

mic EOULInNgG COMMUNnNLEY Structure on the control 
surfaces after six months immersion was dominated either by 
bryozoans or the hydroid Obelia spp. depending on the plate 
examined. Two species of solitary tunicates, Ascidia 
ceretodes and Styela truncata were also censused for the 
first time during the experiment. 

The initial percent cover estimates for the control 
surfaces ranged from 16% to 54%. The mean percent cover for 
the control surfaces was 37% with a standard deviation of 
183. 

The antifouling coated surfaces showed fouling on 
several of the plates with the barnacles Megabalanus 
californicus and Balanus crenatus appearing for the first 
time on these surfaces. Obelia spp. remained the dominant 
fouler on the antifouling coated surfaces. 

THe Similarity values for the control surfaces 
ranged from .48 to .88 with a mean of .70 (Figure 20A). The 
mean Similarity of the antifouling coated surfaces was .99 


(Figure 20B). 
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Figure 20. Similarity Graphs for (A) Non-treated Control 
Surfaces and (B) Anti-fouling Coated Surfaces for Month 6. 
Dashed Lines Indicate Mean Similarity Values. 


64 





2. Computer Simulations 


The computer simulated value of the mean percent 
cover was 33%. The upper 95% confidence limit about the 
mean showed significant decrease out to approximately 20 
plates per group simulated (Figure 21). Note again that 
there was a Significant decrease in the width of the 95% 
confidence limit about the mean as the number of plates per 
group simulated went from ten to twenty. 

The rather high mean similarity value for the sixth 
month control plates and by inferrence the lessened 
interplate variability permitted the 95% confidence limit 
of the fouling cover of 90% of the plates (Figure 22A) to 
be less than that of the fourth month simulations even 
though the mean percent coverage value for the sixth month 
group was higher. The 95% confidence limit of the fouling 
cover of 903 of the plates value decreased until 
approximately the twenty plates per group simulation point 
and then remained fairly constant. 

The standard deviation about the mean percent cover 
that these simulated plates displayed decreased only 
slightly past the 20 plates per group simulated point 
(Figure 22B). Note however that there was a 5% decrease in 
this value between the 10 plates per group simulated point 


and the 20 plates per group simulated. 
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Figure 21. Computer Simulations Using Data From Month 6 
Showing the Expected Value of the Mean Percent Fouling 
Cover (the Mean of the 200 Individual Group Simulation 
Percent Fouling Covers) as a Solid Line and the 95% 
Quantile (Dashed) of the Expected Mean Percent Fouling 
Cover as a Function of Increasing Sample Size. 
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Figure 22. Computer Simulations Using Data From Month 6 
Showing : (A) the Percent Fouling Coverage That One Can Say 
90% of the Samples Will Have Fouling Coverage Less Than or 
Equal to ( With 95% Confidence); and (B) the Standard 
Deviation of Percent Coverage (About the Mean) That These 
Samples Will Have As a Function of Increasing Sample Size. 


67 








G. RESULTS FROM MONTH 7 
1. Experimental Data 

The fouling community structure of the seventh 
month control plates remained dominated by bryozoans with 
the hydroid Obelia spp. still in an important but secondary 
role. The spirorbid worm Circeis armoricana was dominant on 
one of the plates and was common on all the control 
Surfaces. 

The primary fouling organisms on the antifouling 
coated surfaces were the barnacle Balanus crenatus , the 
bryozoan Celloporaria brunnea ,and the hydroid Obelia spp. 
Several of these plates were also fouled by an unknown tube 
dwelling amphipod, possibly of the genus Ampithoe. 

The percent coverage estimates for the control 
Surfaces ranged from 31% to 62% with a mean of 49%. The 
standard deviation of these values was 11.5%. 

The similarity indices between the control surfaces 
ranged from .22 to .88. The mean Similarity value for the 
control plates was .55 (Figure 23A). 

For the antifouling coated surfaces, the range of 
Similarity values was from .8l1 to 1.0. The mean value was 


2N(raqurce 238). 


2. Computer Simulations 


The final iterated value for the simulated mean 


percent fouling coverage was 49%. The upper 95% confidence 
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Pagdre 23. sumitarity Graphs for (A) Non-treated Control 
Surfaces and (B) Anti-fouling Coated Surfaces for Month 7. 
Dashed Lines Indicate Mean Similarity Values. 
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limit about this mean remained fairly constant past the 20 
plates per group Simulated point (Figure 24). 

The Simulated percent cover that one could say 90% 
of the plates would have fouling coverage less than or 
equal to (Figure 25A), decreased slightly until 
approximately the 20 plates per group simulation point and 
then remained fairly constant at approximately 65% 
coverage. The standard deviation that these simulated 
plates had about the mean percent fouling cover showed 
little decrease past the 30 plates per group simulation 


morme (Figure 258). 


H. RESULTS FROM MONTH 8 
i. Experimental Data 

Ene fouling COommunLEeyY Structure on the control 
Surfaces after eight months immersion was dominated nearly 
eeemusavely by bryozoans. The contribution of hydroids to 
the percentage of cover had been reduced on three of the 
surfaces to less than 5%. Serpulid worms, solitary 
tunicates, and barnacles were beginning to emerge as 
important foulers as their increaSingly large size 
prohibited their overgrowth by encrusting species. 

The fouling organisms present on the antifouling 
coated surfaces were primarily the hydroid Obelia spp. and 
barnacles. The spirorbid worm Circeis armoricana was also 


present on several of the plates. 
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Figure 24. Computer Simulations Using Data From Month 7 
Showing the Expected Value of the Mean Percent Fouling 
Cover (the Mean of the 200 Individual Group Simulation 
Percent Fouling Covers) as a Solid Line and the 953% 
Quantile (Dashed) of the Expected Mean Percent Fouling 
Cover as a Function of Increasing Sample Size. 
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Figure 25. Computer Simulations Using Data From Month a 
Showing : (A) the Percent Fouling Coverage That One Can Say 
90% of the Samples Will Have Fouling Coverage Less Than or 
Equal to ( With 95% Confidence); and (B) the Standard 
Deviation of Percent Coverage (About the Mean) That These 
samples Will Have As a Function of Increasing Sample Size. 
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The estimates for the percentage of fouling cover 
ranged from 48% to 97% for the control surfaces. The mean 
of these values was 76% and the standard deviation was 193. 

For the control surfaces, the range of similarity 
values was from .22 to .67. The mean value was .44 (Figure 
26A). 

The similarity values for the antifouling coated 
Surfaces ranged from .75 to 1.0. The mean value was .94 
(Figure 26B). 

2. Computer Simulations 

The computer iterated estimate for the mean percent 
mpoulbing Coverage for the control surfaces was 81%. Note 
Seat the 95¢ coOntidence limit about this value remained 
fairly constant past approximately the 25 plates per group 
Simulation point (Figure 27). 

The high estimate for the simulated mean percent 
cover and the variability permitted by the model forced 
the upper 95% confidence limit about the percent cover that 
90% of the plates would have to remain greater than 96% for 
the entire range of the simulations (Figure 228A). The 
standard deviation about the mean percent cover that these 
plates would have showed a precipitous decrease out to the 
20 plates per group simulated point and then decreased only 
Zamoucmeomtne SO plates per group Simulation point (Figure 


BOO. 
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Figure 26. Similarity Graphs for (A) Non-treated Control 
Surfaces and (B) Anti-fouling Coated Surfaces for Month 8. 
Dashed Lines Indicate Mean Similarity Values. 
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Figure 27. Computer Simulations Using Data From Month 8 
Showing the Expected Value of the Mean Percent Fouling 
Cover {the Mean of the 200 Individual Group Simulation 
Percent Fouling Covers) as a Solid Line and the 95% 
Quantile (Dashed) of the Expected Mean Percent Fouling 
Cover as a Function of Increasing Sample SizZe. 
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Figure 28. Computer Simulations Using Data From Month 8 
Showing : (A) the Percent Fouling Coverage That One Can Say 
90% of the Samples Will Have Fouling Coverage Less Than or 
Equal to ({( With 95% Confidence); and (B) the Standard 
Deviation of Percent Coverage (About the Mean) That These 
Samples Will Have As a Function of Increasing Sample Size. 
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I. RESULTS FROM MONTH 9 
l. Experimental Data 

The bryozoan Watersipora cuculatta dominated the 
fouling assemblages on three of the four control surfaces 
while the hydroid Obelia spp. dominated on the fourth. The 
total number of species represented in the census began to 
Stabilize as the more successful species excluded others in 
the competition for the rapidly diminishing space. 

The dominant fouling organism on the antifouling 
coated surfaces was the bryozoan, Membranipora membranacea. 
This organism is usually found nearly exclusively on the 
Fronds of the giant kelp Macrocystis and it is not known 
maat Chemicak or Yoehner®stimulus attracted it to the 
antifouling coated surfaces. 

The percent coverage estimates for the control 
surfaces ranged from 35% to 92%. The mean percent cover was 
67% and the standard deviation was 25%. 

The similarity values for the plates ranged from 
wEpmEOu.7 4) fOr the control surfaces and from .60 to 1.0 for 
the antifouling coated surfaces. The mean similarity value 
for the control surfaces was .41 (Figure 29A) and was .89 
for the antifouling coated surfaces (Figure 29B). 

2. Computer Simulations 
The final estimate for the mean percent fouling 


cover for the computer simulated plates was 73%. The upper 
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Mmegwire 29... Goimilarity Graphs for (A) Non-treated Control 
Surfaces and (B) Anti-fouling Coated Surfaces for Month 9. 
Dashed Lines Indicate Mean Similarity Values. 
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95% confidence limit about this mean showed little decrease 
past the 20 plates per group simulation point (Figure 30). 
The 95% confidence limit on what one can say 90% of 
the simulated plates will have fouling coverage less than 
Or equal to remained above 95% coverage for the entire 
range of the simulations (Figure 31A). The standard 
deviation that these plates had showed little change past 


the 20 plates per group simulation point (Figure 31B). 


J. RESULTS FROM MONTH 10 
1. Experimental Data 

Bryozoans dominated the fouling assemblages on all 
of the control plates with Watersipora cuculatta occupying 
60% of the space on one of the surfaces. Eleven species of 
bryozoans were identified during the census. 

The similarity values for the plates ranged from 
mOste ./9 £or the control surfaces and £Erom .6/7 to 1.0 for 
the antifouling coated surfaces. The mean similarity value 
for the control surfaces was .47 (Figure 32A) and for the 
antifouling surfaces was .90 (Figure 32B). 

2. Computer Simulations 

The final iterated value for the computer 
Simulation of the mean percent cover was 70%. The upper 95% 
confidence limit about this mean showed little decrease 


Dast the 20 plates per group simulation point (Figure 33). 
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Figure 30. Computer Simulations Using Data From Month 9 
Showing the Expected Value of the Mean Percent Fouling 
Cover (the Mean of the 200 Individual Group Simulation 
Percent Fouling Covers) as a Solid Line and the 95% 
Quantile (Dashed) of the Expected Mean Percent Fouling 
Cover as a Function of Increasing Sample Size. 
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Figure 31. Computer Simulations Using Data From Month 9 
Showing : (A) the Percent Fouling Coverage That One Can Say 
90% of the Samples Will Have Fouling Coverage Less Than or 
Equal to ( With 95% Confidence); and (B) the Standard 
Deviation of Percent Coverage (About the Mean) That These 
Samples Will Have As a Function of Increasing Sample Size. 
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Figure 32. Similarity Graphs for (A) Non-treated Control 
Surfaces and (B) Anti-fouling Coated Surfaces for Month lO. 
Dashed Lines Indicate Mean Similarity Values. 
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Figure 33. Computer Simulatioms Using Data From Month 10 
Showing the Expected Value of the Mean Percent Fouling 
Cover (the Mean of the 200 Individual Group Simulation 
Percent Fouling Covers) as a Solid Line and the 953 
Quantile (Dashed) of the Expected Mean Percent Fouling 
Cover as a Function of Increasing Sample Size. 
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The percent cover that one can say 90% of the 
plates would have coverage less than or equal to (with 953 
confidence) remained above 90% cover for the entire range 
Pewee simulations (Figure 34A). The standard deviation 
these plates would have showed little decrease past the 20 


plates per group simulation point (Figure 34B). 


K. RESULTS FROM MONTH 11 
1. Experimental Data 

The fouling communities on the control surfaces 
remained dominated by bryozoans with Watersipora cuculatta 
dominating on three of the plates and the upright bryozoan 
Bugula californica dominating on the fourth. Of the seven- 
teen species identified on these control surfaces, thirteen 
were bryozoans. 

The predominant fouling organism on those 
antifouling surfaces that showed fouling was the bryozoan 
Membranipora membranacea. The hydroid Obelia spp. and the 
spirorbid worm Circeis armoricana were also in evidence. 

The initial percent coverage estimates for the 
control surfaces ranged from 64% to 98%. The mean value was 
85% and the standard deviation was 15%. 

The similarity values for the plates ranged from 
-28 to .79 for the control surfaces and from .64 to 1.0 for 


the antifouling coated surfaces. The mean similarity value 
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Figure 34. Computer Simulations Using Data From Month 10 
Showing : (A) the Percent Fouling Coverage That One Can Say 
90% of the Samples Will Have Fouling Coverage Less Than or 
Equal to ( With 95% Confidence); and (B) the Standard 
Deviation of Percent Coverage (About the Mean) That These 
Samples Will Have As a Function of Increasing Sample Size. 
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for the control surfaces was .52 (Figure 35A) and for the 
antifouling coated surfaces was .87 (Figure 35B). 
2. Computer Simulations 

The final iterated value for the computer 
Simulations of the mean percent cover was 89%. The upper 
95% confidence limit about this mean decreased by only 7% 
over the entire range of the simulations (Figure 36). 

Quite obviously, with such a large valuefor the 
mean percent cover, the percent cover that 90% of the 
plates will have will be equally large. As can be seen 
(Figure 37A), this value never became less that 97% cover 
over the entire range of the simulations. The standard 
deviation about the mean percent cover that these plates 
would have again showed little decrease past the 20 plates 


per group simulation point (Figure 37B). 


86 








D 
SIMILARITY 
0.50 0.75 1.00 


0.25 


0.00 





D18F 018B O2oF 025B 
PLATE DESIGNATION : 


es 
SIMILARITY 





iemenipencen 205 BohwbSb Bi7k Bl/B 02h C2B C23F Cis 
PLATE DESIGNATION 


Figure 35. Similarity Graphs for (A) Non-treated Control 
Surfaces and (B) Anti-fouling Coated Surfaces for Month ll. 
Dashed Lines Indicate Mean Similarity Values. 
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Figure 36. Computer Simulatioms Using Data From Month 11 
Showing the Expected Value of the Mean Percent Fouling 
Cover (the Mean of the 200 Individual Group Simulation 
Percent Fouling Covers) as a Solid Line and the 953 
Quantile (Dashed) of the Expected Mean Percent Fouling 
Cover aS a Function of Increasing Sample Size. 
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Figure 37. Computer Simulations Using Data From Month 11 
Showing : (A) the Percent Fouling Coverage That One Can Say 
90% of the Samples Will Have Fouling Coverage Less Than or 
Equal to ( With 95% Confidence); and (B) the Standard 
Deviation of Percent Coverage (About the Mean) That These 
Samples Will Have As a Function of Increasing Sample Size. 
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V. CONCLUSIONS AND RECOMMENDATIONS 


A. DISCUSSION 

The choice of an appropriate sampling strategy for the 
study of fouling organisms in Monterey Bay has been shown 
to be dependent upon the type of information desired. 

By analyzing the width of the upper 95% confidence 
limit about the expected value of the mean percent fouling 
cover, estimates were made concerning the optimum number of 
plates to be deployed. This optimum number was found by 
determining that point where the addition of additional 


Plates Mad a negligible effect upon the width of the 


confidence interval. In nearly all the cases simulated, the 
value for the optimum number of plates to deploy appeared 
to be approximately twenty. 

A similar procedure was followed to estimate the 
optimum number of plates to deploy to minimize the standard 
deviation about the mean percent fouling cover that 90% of 
the plates would have. The optimum number of plates to 
deploy to satisfy this requirement was again estimated to 
be approximately twenty for most of the simulations. 
However, the results from several cases suggested that 


thirty plates would be a more appropriate sample size. 
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pacedmonmeenese results, 1t is concluded that twenty 
plates is probably the minimum number of plates that should 
be deployed to obtain accurate estimates of the mean 
Pemoent fouling cover of a group of plates for this 
locality. Thirty plates per group is probably a more 
appropriate number of plates to deploy to insure that 
accurate results are obtained for groups that display 
heightened variability in the individual plate fouling 
coverage estimates. 

The experimental results showing the negligible effect 
of the addition of more plates on what one could say 90% of 
the plates would have fouling coverage less than or equal 
to (with 95% confidence) lead to the conclusion that the 
mime renmt Vartability of Fouling populations is 
Significantly greater than previous studies indicated. This 
means that while twenty to thirty plates are probably 
sufficient to resolve ambiguities concerning the mean per- 
cent fouling cover, this number is clearly insufficient to 
ascertain with any high degree of confidence the amount of 
variability of a large segment of the population. For 
example, if it is desired to ascertain with 95% confidence 
the variability of 90% of the fouling population, it must 
be understood that this will require the committment of 


Substantial resources to the study. 
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In addition to the development of an appropriate 
fouling community sampling strategy for Monterey Bay, moce 
far reaching conclusions can be drawn concerning the 
applicability of the procedures used in this thesis to 
other locations. Since the computer modelling procedure 
wsead Mm this study dealt with the extension of 
the experimentally obServed variability and contained no 
Site-specific parameters, it is believed that the procedure 
can be directly applied to the study of fouling community 
variability at any desired geographical location. This 
means that the computer programs developed for this study 
can be coupled with archived fouling coverage data from 
amectce, GCepth, season, and so on to provide information 
on how best to sample the fouling populations. 

The final conclusion drawn from this study is that the 
bootstrap method of computer intensive statistical analysis 
has profound implications in the study of other biological 
problems. It is believed that this method could be applied 
to a wide range of other data intensive biological areas 
including fisheries management, larval settlement Studies, 
the pelagic distributions of plankton and nekton, and 


growth rate studies to name just a few. 


B. RECOMMENDATIONS FOR FURTHER RESEARCH 
The model developed in this thesis should be applied to 


a number of other geographical locations to determine the 
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differences in fouling variability in terms of latitude, 
longitude, depth, environmental stresses, and a host of 
Sener biological forcing functions. Since the computer 
program already exits and the only required inputs are 
archived fouling coverage data, this should present no 
insurmountable difficulties. 

Once such studies have been completed, the development 
of empirically derived sampling strategies for any location 
could be attempted. 

Finally, the use of the bootstrap or other computer 
intensive statistical techniques should be vigorously 
investigated in other areas of biological interest. It is 
believed that the use of these techniques might well 
provide answers to a wide range of biological problems that 


have so far proved intractable. 
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APPENDIX A 


MICRON 22 ORGANO-METALLIC 
POLYMER ANTIFOULING PAINT 


INGREDIENTS PERCENT BY WEIGHT 
Active: 
BES Gem OUty eein) Oxide liven, 
Cuprous Thiocyanate Tacs 
Inert: ees 
100 


Elemental Tin 4.4% 
Elemental Copper 8.9% 


Potmemcontains 1.1 lbs of Bis({trfibutyltin) oxide per 
gallon and 1.6 lbs of Cuprous Thiocyanate per gallon. 
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APPENDIX B 


NAVY STANDARD FORMULA 121 
RED VINYL ANTIFOULING PAINT 


INGREDIENTS PERCENT BY WEIGHT 
Cuprous Thiocyanate 7 Ones 
Rosin On > 
Vinyl Resin ae) 
Tricreysl Phosphate 2.4 
Methyl Isobutyl Ketone Creel: 
Xylene 56 
Antisettling Agent 4 


Source: Department of the Navy Specification MIL-P-1L5931C, 


Paint, Antifouling,Vinyl(Formula Numbers 121 and 129) 
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APPENDIX C 


NAVY STANDARD FORMULA 170 BLACK 
CAMOFLAGE ANTIFOULING PAINT 


INGREDIENTS PERCENT BY WEIGHT 
Vinyl Resin io 
Els trloultyLetin) Oxide 4.2 
TELOUEyE Cin BLueride bol 
Carbon Black 251) 
Titanium Dioxide =o 
Ethylene Glycolmonoethyl B40 
Ether Acetate 
Normal Propanol ees li 
Normal Butyl Acetate 43.2 


Source: Department of the Navy Military Specification DOD- 


P-24588,Paint,Antifouling,Vinyl,Camoflage(Formula numbers 
Wine pki ana 173), 2 May 1979. 
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APPENDIX D 
ZYNOLYTE EPOXY RUST MATE PAINT 
INGREDIENTS PERCENT BY WEIGHT 
Non-Volatile (58.4% of total) 
Pigments 43.4 
Vehicle 56.6 


Epoxy and Menhaden 
Alkyd Resins 


HOOT 

Volatile (41.6% of total) 
Exempt Mineral Spirits 98.0 
Aromatic Hydrocarbons 250 
n8O) 0: 


Source: Product ingredient breakdown on label 
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APPENDIX §E 
A DISCUSSION OF THE METHOD OF MAXIMUM LINELIHOOD AND ITS 
INCORPORATION INTO A MODEL FOR FOULING COVER 


A. GENERAL 

The following model was used to describe the properties 
of the proportion of a plate that was fouled. Each plate 
(t) has a random proportion of fouling,P,. Given Pi the 
number of the 100 censused sites that are fouled 1 (Si), on 
the tth plate has a binomial distribution: 
100-k oe 


OL eee) 


P{S, = Jie « _ a - 


for k = 0,1,2.,100 independent of the other plates. The 


proportion of fouling,P,, is assumed to have the form 


€ 
a5 (2) 


* ase’ 
where é, has a normal distribution with mean yp and 
variance oa” independent of the other plates censused at 
the same time. The mean wu and variance o* will in general 
be a function of the amount of time the plate is submerged 
though the many variables involved in population dynamics 
will keep the function from being linear. 


Assume N plates are inspected at a time with the result 


ehlat S, of the censused sites are fouled on the ttf plate. 
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The likelihood function for the model is: 





N 
S 100-S 
2 100 c ae 
[as | ote, 3 wee) % ) P,. (jl - eee (3) 
ter: 
where , 
ee 
at 2 
2¢ ) as 
fe, 3,0°) = = O for -0 < e, < (4) 


ov2T 


Note that L is the probability of observing Sp oreer 
fouled sites on each of the plates. 
The method of maximum likelihood is to find those 


values of UM ,oa, and € fOree = 72,...N Which Maximize 


t v 
L; that is, those values which maximize the probability of 


observing the outcome. These values also maximize the log 


likelihood function, in L, where 


N P 

Ce, - us eG 

Ys inl =o fe eae + S,anP, - tno + (100-S,)en(1-P.)} (5) 
al 


[+ constant ] 


To find the maximum of a the partial derivatives of ig 
with respect to yw, o and €& are set equal to zero. Tides 


results in the equations: 
N 


ae . ae 
du 3° 
t=1 


solving for mu results in: 
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Similarly N 








2 
ae a ee 
cl ‘oe >) a 
0 e 
t=l / 
solving for o* results in 
N 
(7) 
a2 _1 : : 2 
OFF Ce, - Ul) 
iy 
and ( ¢ yp 
a ie 
O t t a ik: iG 
ae. K bee Ss] 2S OES) [-+~—-] =0 (8) 
de, o i Ee, nc 1 GERD) Oe, 
where oP, : te ae 
oe, - E, [ E - Eo ' PCP, 
lte (lte ~) 
Equation (8) simplifies to 
€ 
te 
Z 9 
ei ‘Cue6) a =— 40 — 5,0 = 0 a 
1E . 
ite 


B. SOLUTION OF THE LOG LIKELIHOOD EQUATION 


A recurSive method was used to solve the system of 
equations (6),(7), and (9). For each month, the 100 sites 
were censused on each of the four untreated plates and S&S, 


the number of fouled sites was determined. Initial values 


EOE PB, »€, us and a" were determined as follows 


“A 


Each) = S,/100 


ae (=PNOT) 


100 





ae = P _p =EPSN 
€ (0) = antP,(O)/C1 P,.(0))] ( SNOT) 
u 
ap 5 ) € (0) (=MUNOT) 
ca 
u 
A 7 1 : ) ~ i 2 
ea 


The values ino} and Go }were then used in equation (9) and 
the equation solved for cae for t = 1,2,3,4. Since 
equation (9) is transcental, Newton's method was used to 
solve for €(1). Newton's method iS an iterative procedure 


in which the (n+1)St member of the iterative sequence (Xp) 


a £CX_) 


where in our case 
x 
n a a a9 
£(X.) =X, + Sy (100 57(0)) = fi(0) - S,0°(0) 
n 2 


x 
n 





lt+e 
and 


x x > 
eee ee 2/Gre)) ialog cca) 
hi 


Equation (10) was iterated until |x ay | <1lx 107° 
Masia n 
~» 
The é (1), t = 1,2,3,4 were then used to compute iterated 
A A. A 
values of wU, OO, and a as follows: 


A ¢.(1) a) 
P.(1) =e Se arco sn) 


OL 





u 


A pee A 
yl) =a >: e,(1) 


t=1 
un 
é(1) = 2 ys (2,() - wD)? 
t=1 
and equation (9) was then solved with i = i (1) and G 


=o (1) for E(2), t = 1,2,3,4. The iterative procedure 
was continued until the achievement of a tolerance value of 
—6 \ » in 

(1x10 ) calculated by |u(k) - u(k+l)| indicated the system of 
equations had converged. The final iterated values for 
a a A2 A. 

P. (=PNEW), €. (=EPSNEW), o° (=SIGNEW), and UM (=MUNEW) 
were the values that maximized the likelihood function for 
the model. SIGNEW (@) and MUNEW (1) were then used as 


initial input values in the bootstrap simulations of the 


experiment (Appendix F). 
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APPENDIX F 
BOOTSTRAP COMPUTER SIMULATIONS 
A. AN EXPLANATION OF THE BOOTSTRAP METHOD OF COMPUTER 
SIMULATIONS OF RANDOM PROCESSES AND ITS INCORPORATION 
INTO AN ASSUMED STOCHASTIC MODEL FOR FOULING COVERAGE 
ie Discussion 
The common statistical tools utilized in the study 
Semevotrculamg have as their basis the simplifying 
assumption that the data collected from the analysis of 
such communities can be described by a normal or Gaussian 
Guscretoution. Ehat 1S to say, it is assumed that 
Fluctuations in the values of some experimentally observed 
parameter are scattered symmetrically about the true value 
of the parameter. It is further assumed that the larger the 
difference between the the true value and the observed 
value of the parameter, the less likely it is that the 
value will be observed experimentally (Diaconis and 
Efron,1983). Many years of experience using these 
assumptions have shown that even if the dataare only 
approximately or "pseudo" normal, the Gaussian theory still 
works quite well. If however, the data do not satisfy the 
requirements for the assumption of normality or,if the 
sample size is such that the various tests used to check 


for normality can only give ambiguous results, it is clear 
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tmiae tie results of statistical techniques based on the 
assumption of normality will be unreliable. 

Recent developments in the use of computer- 
intensive techniques for statistical analysis, particularly 
the invention of the bootstrap technique (Efron,1977), have 
enabled the computation of various statistical parameters 
without the necessity of assuming a Gaussian distribution. 
This technique has also enabled the computation of those 
Statistics which do not have a simple analytical formula. 
Prior to the advent of large main frame computers, the 
Q@z£Eiculty of finding numerical solutions to non-linear 
problems forced statistical methods to concentrate on those 
Statistical models and procedures for which analytical 
results could be obtained. These models and procedures did 
give useful large sample size results for the common 
Statistics such as the mean and variance of a population. 
However, they ignored other important statistical questions 
that did not have analytical formulas; such as, the degree 
of variability in an estimate due to the sample being of 
finite size or the confidence limits about estimates froma 
finite sample (Diaconis and Efron,1983). In general terms, 
the bootstrap method consists of coupling a probabilistic 
model with the data gathered experimentally from one sample 
of size N to generate a large number of simulated samples 


of size N. These simulated samples are then analyzed to 
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determine the variability of estimates of the true values 
of the statistics of the population. 
Zee SOCeadune 

The model described in Appendix E was used in a 
Simulation study of the variability of the estimate of the 
Fraction of plate coverage. The model was used to simulate 
200 groups of fixed numbers of plates. The inputs into the 
Simulation were the values for the maximum likelihood 
estimates of the mean (MUNEW i ) and variance 


(SIGNEW oa POr the monthlyaeepsilon me values 


+) 
determined from the analysis of the fouling cover on the 
four untreated plates. 

The stochastic model used in this casSe assumed a 
normal distribution for the epsilon (€,) values ( that is 


A 


€é,~ N(i,6)). Using available computer software, the 
required number of normally distributed random numbers with 
mean equal to zero and variance equal to unity were 


generated for 200 groups of the following number of plates: 


Number of Simulated Required Number (TI) 
Plates per Group of Random Numbers 
2 I = 400 
a I = 800 
5 I = 1000 
10 I = 2000 
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Number of Simulated Required Number (TI) 


Plates per Group of Random Numbers 
15 I = 3000 
20 I = 4000 
30 I = 6000 
40 I = 8000 
80 I =16000 


Eacn of the (I) computer generated random numbers 
from the standard normal distribution were then transformed 
into random numbers(n,;) with a normal distribution with 
mean i and standard deviation @ by multiplying each 
ithyariable by @ (the Square root of @G) and adding i : 
Soe Fandom Number n; was transformed to give PINIT(I) 
defined as: 

PINIT(I) = eNY(1 + ee) , the average proportion of 
the ith plate fouled. 

The initial simulated percent cover for the jth 
plate was determined by calculating a random number having 
the properties of a binomial distribution with probability 
(P) equal to PINIT(I) and N equal to 100. The resulting 
variable was termed PNOT(I). The variables EPSNOT(I),MUNOT, 
SIGNOT, PNEW(I),EPSNEW(I), MUNEW, and SIGNEW were then 


calculated using the method of maximum likelihood described 


in Appendix E, 


106 








MONTH 


MONTH 


MONTH 


MONTH 


MONTH 


MONTH 


APPENDIX G 


TABULATED MONTHLY PERCENT FOULING COVERAGE VALUES 
FOR THE NON-TOXIC CONTROL SURFACES 


PLATE 


PLATE 
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mJ NF 
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eee | 
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PLATE £ 
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| eae 
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FRACTIONAL COVERAGE 
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APPENDIX 4H 


Piette sohoSTLE SPECIES IDENTIFIED BY THE RANDOM POINT 
CENSUS AND THE MONTHS THEY WERE PRESENT ON THE NON-TOXIC 
CONTROL SURFACES 


MONTH NUMBER 
ge) (Ora tin 6 


ORGANISMS 2 


J 
|\o 
_ 
© 
a) 


Protozoa: 
Folliculina spp. ho xk UX x ao X 
x 


Ephelota gemmipara 


Coelentrata: 
Obelia Spp. Ne eh As om aX 


Hydractinia spp. X 
Anthopleura spp. X 


Ectoprocta: (Bryozoans) 
Bugula neretina 
Bucula (Gali fornica 


Watersipora cucullata 
Hippothoa hyalina 
Celloporaria brunnea 
Cryptosula pallianasa 
Schizoporella unicornis 
Muecreporella ciliata 
Microporella californica 
Membranipora membranacea 


Membranipora serilamella X 
Unknown bryozoan #1 


Unknown bryozoan #2 X 


PS PS 


vs 
vS 
rs 


PS OS OS OS 
MS OS OS OS 
PS OS PS OS PS OS OS OS 
PS OS PS OS XS OS 
PS OS OS OS OS 
VS OS PS 
PS OS OS OSS OS OS OM 
PS PS OS OM PS 
POS MS OM OM 
PS PS PS OS PS 


PSS PS OS 
VS PS 


Annelida: 
Circeis armoricana <a x 


Janua nipponica 
Pileolaria potswaldi 


Protolaeospira exima Xx 
Serpula vermicularis x aX ms 


Anatides groenlandic groenlandica Xx X 


Arthropoda: 
Balanus crenatus Xx X 
Megabalanus californicus x X xX X xX 
Amphipod (unknown) x 





PSS XS 
PS PS PS 


Mollusca: 
Mytilus edulis X 


eS 








ORGANISMS 


Echinodermata: 
strogylocentrotus spp. 


Chordata: 
Ascidia ceretodes 
ptyela truncata 
Pyura haustor 


La) 


MONTH NUMBER 


Oe ee 
X one 
X 


10 








LIST OF REFERENCES 


Boyd,M.,1972. Fouling Community Structure and Development 
in Bodega Harbor, California. Doctoral Dissertation. Unive- 
rsity of California,Davis. 


Diaconis,P., and Efron,B.,1983. Computer-Intensive Methods 
in Statistics.Scientific American.,v.248,No.5:116-130. 


Eiron,5.,1979. Computers and the Theory of Statistics: 
Thinking the Unthinkable.In SIAM Review,v.21,No.4.0cto- 
ber,1979:460-480. 


Fisher,E.C., Birbaum,L.S.,Depalma,J.,Muracka,J.S.,Dear H., 
and Wood,F.G. 1975. Survey Report: Navy Biological Fouling 


and Biodeterioration. Naval Undersea Center Report NUC-TP- 
456. 


Fraser,C.M.,193/7.Hydroids of the Pacific Coast of Canada 
and the United States, University of Toronto Press.297pp. 


Haderlie,E.C. 1974. Growth Rates, Depth Preference, and 
Ecological Succession of Some Sessile Marine Invertebrates 
in Monterey Harbor. Veliger.v.17 (supplement) :1-35. 


Kelley,P.R.,1981. Scanning Electron Microscope Observations 
of Marine Microorganisms on Surfaces Coated with 
Antifouling Paint. Master's thesis, Naval Postgraduate 
Senoolemonterey, California, U.S.A. 


Knight-Jones,P.,and Knight-Jones,E.W.1979. Spirorbidae 
(Polychaeta Sedentaria) from Alaska to Panama. J. Zool. 
Lond.v.189:419-458. 


Mook,D.,1976. Studies of Fouling Invertebrates in the 
PictanimiumyetemSilLletin of Marine Science.v.26:610-615. 


Morris,R.H.,Abbott,D.P.,and Haderlie,E.C.,1980. Intertidal 
imiwcmuciomeaecs Of California.-Stanford University 
Press. Pp. 


Osburn,R.C.,1952.Allen Hancock Pacific Expeditions v.14, 
parts 1,2,and 3 , Univ. of Southern California Press.841pp. 


Osman,R.W.,1977. The Establishment and Development of a 


Marine Epifaunal Community. Ecological Monographs.v. 
A7:pp 37-63. 


ileal 








Schoener,A.,Long,E.R.,Depalma,J.R.,1978. Geographic 
Gwar brawvon im ArELi£Eicial Island Colonization Curves. 
Ecology.v.47,No.2:36/7-382. 


Schoener,A., and Greene,C.H.,1980. Variability Among 
Identical Fouling Panels in Puget Sound, Washington, 
U.S.A.In Proceedings of the Fifth Internationl Congress on 
Marine Fouling and Corrosion.Barcelona ,Spain.May 1979. 


Sm th, R.I.,and ear oton,o. la,1975. Eds. Light's Ss Manual, 


3rd ed.,University of California Press. 716pp. 


Sutherland,J.S.,1974. Multiple Stable Points in Natural 
communities. American Naturalist.v. 108:859-873. 


Whitaker,R.H.,1952. A Study of Summer Foliage Insect 
Communities in the Great Smoky Mountains.Ecological 
Monographs.v.22,No.l:pp 1-42. 


Woods Hole,1952. Marine Fouling and its Prevention. U.S. 
Naval Institute, Annapolis, Md. 


eZ 











Peete DiStTRIBUTION LIST 


No. Copies 


Defense Technical Information Center Z 
Cameron Station 
Alexandria, Virginia 22314 


Library, Code 0142 2 
Naval Postgraduate School 
Monterey, California 93940 


Professor Robert J. Renard, Code 63Rd iL 
Department of Meteorology 

Naval Postgraduate School 

Monterey, California 93940 


Professor Christopher N. K. Mooers, Ni 
Code 68Mr 

Department of Oceanography 

Naval Postgraduate School 

Monterey, California 93940 


Professor Eugene C. Haderlie, Code 68Hc l 
Department of Oceanography 

Naval Postgraduate School 

Monterey, California 93940 


Professor Patricia A. Jacobs, Code 55Jc i 
Department of Meteorology 
Naval Postgraduate School 


Professor Donald P. Gaver, Code 55Gv 1 
Department of Operations Research 

Naval Postgraduate School 

Monterey, California 93940 


Loewe pavia &. Martin 1 
Route 3, Box 667 
Yakima, Washington 25314 


Director 1 
Naval Oceanography Division 

Naval Observatory 

34th and Massachusettes Avenue NW 

WashuncgvonneD. Cc. 20390 


eS 


roe 


vee 


Ve 


or 


14. 


eye 


16. 


ye. 


ae. 


or 


Commander 

Naval Oceanography Command Central 
NSTL Station 

Bay St. Louis, Missouri 39522 


Commanding Officer 

Naval Oceanographic Office 
NSTL Station 

Bay St. Louis, Missouri 39522 


Commanding Officer 
Fleet Numerical Oceanography Center 
Monterey, California 93940 


Commanding Officer 

Naval Ocean Research and Development 
Activity 

NSTL Station 

Bay St. Louis, Missouri 39522 


Commanding Officer 

Naval Environmental Prediction Research 
Facility 

Monterey, California 93940 


Chairman, Oceanography Department 
U.S. Naval Academy 
Annapolis, Maryland 21402 


Chief of Naval Research 
800 N. Quincy Street 
RMaMInNgton, Virginia 22217 


Office of Naval Research (Code 480) 

Naval Ocean Research and Development 
Rete) ty. 

NSTL Station 

Bay St. Louis, Missouri 39522 


Commander 

Oceanographic Systems Pacific 
Box 1390 

Pearl Harbor, Hawaii 96860 


Mrs. Anne Harrington 


Hopkins Marine Station 
Pacific Grove, California 93950 


114 








Zur 


Ze 


Dr. Amy Schoener 
Department of Oceanography 
University of Washington 
Seattle, Washington 98195 


Mrs. Jean Montemarino, Code 2844 

David W. Taylor Naval Ship Research 
and Development Center 

Annapolis, Maryland 21402 


> 














201646 





201646 


Thesis 
M3573 Martin 
Ca The use of computer 


intensive statistical 
modeling in estimating 
the variability of 
marine fouling communi- 
ties. 





